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I.  INTRODUCTION 


It  is  well  known  that  a  slender,  high-performance  ballistic  reentry 
vehicle  is  susceptible  to  roll  resonance  caused  by  small  asymmetries  in  mass 
and  configuration.^-^  The  vehicle  can  exhibit  large  roll-rate  excursions  be¬ 
cause  of  its  low  roll  moment  of  inertia  and  the  extreme  aerodynamic  pressures 
during  atmospheric  entry.  The  vehicle  Is  most  susceptible  to  a  body-fixed 
trim  asymmetry  with  an  orthogonal  component  of  radial  c.g.  offset.  Lift 
caused  by  the  trim  angle  of  attack  acts  on  the  c.g.  offset  moment  arm  to 
produce  a  roll  torque  that  can  spin  the  vehicle  rapidly  into  resonance.  A 
trim-angle-of-attack  asymmetry  on  the  order  of  I  deg  or  less  in  conjunction 
with  a  radial  c.g.  offset  on  the  order  of  tens  of  thousandths  of  an  inch  can 
cause  large  roll  rate  excursions  sufficient  to  spin  the  vehicle  into 
resonance.  The  trim  angle  of  attack  can  be  amplified  in  resonance  by  a  factor 
of  10  to  15  or  more,  depending  on  the  altitude  at  which  resonance  is 
encountered.  As  a  preventive  measure,  vehicles  must  be  accurately  balanced  to 
minimize  radial  displacement  of  the  center  of  gravity  from  the  aerodynamic 
axis  of  symmetry.  In  some  cases,  roll  control  is  also  required. 

Described  here  Is  a  method  of  utilizing  roll  resonance  in  a  controlled 

manner  for  drag  modulation.  Large  angle-of-attack-induced  drag  can  be  used  to 

decelerate  a  test  vehicle  for  recovery,^  or  drag  modulation  can  be  used  to 

£ 

compensate  for  drag  uncertainty  in  order  to  control  range  errors.  The 
control  system  consists  of  a  moving-mass  roll  control,  strapdown  motion 
sensors,  and  a  built-in  trim-angle-of-attack  asymmetry.  During  roll  reso¬ 
nance,  the  roll  rate  is  approximately  equal  to  the  undamped  natural  pitch 
frequency,  which  can  become  quite  large  for  a  ballistic  reentry  vehicle  near 
peak  dynamic  pressure.  However,  the  frequency  at  which  the  roll  rate  oscil¬ 
lates  about  the  critical  frequency,  which  determines  the  response  frequency  of 
the  moving-mass  control  system,  is  considerably  lower  than  the  pitch  fre¬ 
quency.  This  feature  of  the  control  system  as  well  as  the  relatively  small 
mass  required  for  c.g.  control  greatly  simplifies  the  practical  implementation 
of  the  system.  A  feedback  law  is  derived  from  a  linear  approximation  to  the 
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moment  equations  of  motion  near  resonance.  The  control  system  as  applied  to 
recovery  is  demonstrated  with  a  digital-computer  simulation  of  the  nonlinear 
moment  equations  of  motion. 


II.  CONTROL  ANALYSIS 


The  small-angle  equations  of  rotational  motion  in  terms  of  classical 
Euler  angles  or  wind-referenced  coordinates  can  be  written 


•  •  n  •  n  •  r\ 

0  +  (a)  —  )0  +  U 0  =  CO  T  COs(<(>  +  4>  q) 


(1) 


,  •  •  •  •  2 

t-  (4>0 )  +  0\j)  +  v4»0  =  a)  t  sin (4>  +  <J)q) 


(2) 


p  =  -p0c  sin<j> 


(3) 


in  which  t  is  a  body-fixed,  nonrolling  trim  angle  of  attack  oriented  at  an 
angle  <j>Q  with  respect  the  plane  of  a  radial  c.g.  offset  c  (Fig.  1).  The  roll 
angle  <j>  is  a  measure  of  the  orientation  of  the  plane  of  c.g.  offset  with 
respect  to  the  wind  plane.  We  will  consider  the  case  in  which  the  trim  is 
orthogonal  to  the  plane  of  c.g.  offset  (<j>Q  =  90  deg),  which  will  cause  a  rapid 
spinup  into  resonance.  We  assume,  a  priori,  the  quasi-steady  condition 

xJj  «<  ou  =  const  (4) 

which,  with  the  small-angle  definition  of  p 

•  • 

p  =  4>  -H  (5) 


gives,  for  Eqs.  (2)  and  (3), 
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$  +  n8c  sin<j> 
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We  further  assume  that  the  vehicle  spins  through  resonance  so  that  the  body- 
fixed  trim  plane  is  oriented  nominally  180  deg  from  its  nonrolling  orienta¬ 
tion,  $  =*  -90  deg  (Fig.  2),  and  define  a  new  roll  angle  e 


e 


( 


If  we  assume  that  e  is  a  small  angle,  Eqs.  (6)  and  (7)  can  be  written 


U)T 

2 


e 


e  -  n0c  =0  ( 


As  the  vehicle  spins  into  resonance,  the  moving-mass  controller  will  adjust 
the  center  of  gravity  offset  to  cause  a  roll  response  that  will  control  the 
angle  of  attack  or  normal  acceleration  according  to  some  control  law.  We 
select  a  control  law 


K.e  +  K,.(0  -  8  )  +  K,,G 


in  which  0^  is  a  command  value  of  angle  of  attack  (or  normal  acceleration). 
Equation  (11)  substituted  in  Eq.  (10)  gives 


£  +  Kxe  +  1^(9  -  6  )  +  1^0  -  0  (12) 

which,  with  Eq.  (9),  yields  two  coupled  control  equations  in  0  and  e.  We  can 
solve  for  0  from  Eqs.  (9)  and  (12)  by  the  use  of  Laplace  transforms,  which 
gives 

[s3  +  (Kx  +y)s2  +  (JKj  +~2-K3)s 

(13) 

+-K2]0(s)  -^K29c(s) 

Consider  the  response  to  a  step  command  0^(0  “  9c«  Equation  (13)  can  then  be 
written 


0(s) 


K0 

_ c _ 

s(s  +  a)  (s  +  b)  (s  +  c) 


(14) 


where 


(s  +  a)(s  +  b)(s  +  c)  a 

s3  +  (Kx  +  ^-)s2  +  YLl  +  ~  K3)s  +  K 


(15) 


and  K  =  tux / 2 .  We  can  obtain  a,  b,  and  c  from  a  root  locus  in  which  the 
transfer  function  G  is 


s[s2  +  (Kj  +  j)s  +  %  Kx  +  ^  K3] 


Two  possible  root  loci  are  shown  in  Fig.  3,  depending  on  the  coefficients  of  8 
in  Eq.  (16).  For  a  stable  solution,  either  three  negative  real  roots  or  one 
negative  real  and  two  complex  conjugate  roots  are  possible.  The  inverse 
Laplace  transform  of  Kq.  (14)  yields  the  time  response  to  the  step  command  0 


b(b  -  a)  (b 
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III.  APPLICATION  TO  RECOVERY 


The  large-angle  pitch  and  yaw  equations  equivalent  to  Eqs.  (I)  and  (2) 
can  be  written"* 


0  -  ijj  sin0  cos0  +  v0  = 

2  /  i  i  ,  \  CN(6)  qSxst 

0)  T  COS($  +  4>q  ) - | - 


(18) 


*  . 

a  1 

q 


(<ji  sin0)  +  0\j>  cos0  -  yp0  +  ij)  sin9 


-  CN  pp0  »  ui  t  sin( 4>  +  (f>0) 
a 


C^(6)  'J'  sln^  cos0 
a _ 

V9> 


(19) 


in  which  the  lift  coefficient  0^(0)  and  the  normal  force  coefficient  C^(0)  can 
be  approximated  by  the  sharp  cone  Newtonian  relations. ^  For  0  <  o 


C^(8)  =  cos  o  sin20 


C.(0)  =  2  sin^o  +  (l  -  3  sin^o)  sin^0 


(20) 


(21) 


and  for  0  >  a 


CN(6)  «  (cos^a  sin  29) [  ) 


+  (co5_8Kx  +  2)] 


(22) 
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P4rp 


(23) 


C  (0)  -  (iL±J.Zl)  [2  s in2 a  +  (l  -  3  sin2a)] 

A  7T 

3 

+  cos 3  sin20  sin2a 

where  a  ■  cone  half-angle,  X  =  tan  a/ tan  0,3=  sin  ^X,  and 


C  (0)  -  C.(0)  cos0  +  C„(0)  sin©  (24) 

DA  N 

C.  (0)  -  C  (6)  cos0  -  C.(0)  sin0  (25) 

L  N  A 


We  can  estimate  influence  of  the  angle  of  attack  and  drag  on  the  trajectory 
from  Eqs.  (3),  (18)  and  (19)  in  conjunction  with  the  trajectory  equations-* 


du 

dt 


CDqS 


m 


g  siny 


(26) 


dh 

dt 


-u  siny 


(27) 


dy 

dt 


R, 


“ _ ) 

,  +  h' 


cosy 


(28) 


The  moving-mass  controller  is  simulated  by  the  c.g.  offset  relation  given  in 
Eq.  (11)  in  which  £  «  -$  according  to  Eq.  (8).  The  tuncfion  of  the  control 
system  is  to  limit  the  normal  acceleration  to  some  acceptable  value,  while 
providing  sufficient  angle-cf-attack-induced  drag  to  decelerate  the  vehicle 
from  hypersonic  velocity  at  recovery-initiation  altitude  to  a  soft  landing. 
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IV-  NUMERICAL  EXAMPLE 


The  equations  of  motion  were  solved  numerically  for  a  recovery  simulation 
in  which  a  ballistic  reentry  vehicle  is  spun  into  resonance  at  an  altitude  of 
15  kft.  The  moving-mass  controller  is  utilized  to  control  the  growth  of  angle 
of  attack  in  order  to  limit  the  normal  acceleration  to  a  prescribed  maximum 
value.  The  angle  of  attack  is  related  to  normal  acceleration  AN  according  to 


V0) 


mg  An 
qS 


(29) 


>,  hich,  for  small  angles,  can  be  written 


mg  x  . 

o  mg  ,  °  st 

0  “  ~  \  =  2t  ^ 
0)  I 


CN  “S 


a 


(30) 


The  control  law  Eq.  (11),  expressed  in  terms  of  normal  acceleration  and 
normal  acceleration  command  Ajq  is 


.  .  mg  x 

-  nTilii  K*  +  K2 - 27  (AN  -  V  -  K39! 

1  0)  I  c 


(31) 


The  control  configuration  is  such  that  the  plane  of  c.g.  control  is  orthogonal 
to  the  plane  of  trim  asymmetry  (<j>g  -  90  deg  in  Fig.  1).  The  vehicle  having 
the  aerodynamic  and  mass  properties  listed  below  is  subjected  to  a  suddenly 
applied  trim-angle-of-attack  asymmetry  t  =  1.5  deg  with  an  initial  c.g.  offset 
of  0.050  in. 
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hQ  -  15  kft 


Jo  =  15.75  kft/sec 


Am  =  400  3 


d  =*  0.958  ft  “  26  deg 

I  =  3.70  slug-ft2  a  ■  6  deg 

I  ■  0.140  slug-ft2  p  “50  rad/sec 

x  o 

m  =3.79  s  lugs  =  400  g 

c 

S  «  0.721  ft2 
xgt  “  0.468  ft 

A 

As  the  vehicle  spins  into  resonance,  the  c.g.  offset  is  controlled  according 
to  Eq.  (31),  in  which  the  feedback  gains  are  as  follows: 


t  <  1.33  sec 


t  >  1.33  sec 


22469 


The  gains  were  estimated  from  the  linear  results  of  Eqs.  (14)  through  (16)  to 
give  a  stable  solution  and  were  changed  once  to  accommodate  the  large  change 
in  dynamic  pressure  and  vehicle  characteristic  frequencies  as  the  vehicle  de¬ 
celerates.  The  large  angle-of-attack  aerodynamics  are  calculated  from 
the  Newtonian  approximations,  Eqs.  (20)  through  (25),  and  the  pitch  damping 
derivative,  ,  is  altered  as  a  function  of  angle  of  attack,  as  shown  in  Fig. 

4,  to  approximate  the  destabilizing  effects  of  vortex  shedding  at  large 

8  9 

incidence  described  in  NASA  Ames  experiments.  * 
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Results  of  the  computer  simulation  are  shown  in  Figs.  5  through  7.  The 
roll  rate  response  relative  to  the  critical  roll  rate  is  shown  in  Fig.  5a,  the 
c.g.  offset  movement  is  shown  in  Fig.  5b,  and  the  windward-meridian  angle  4>  is 
shown  in  Fig.  5c.  The  normal  acceleration  response  to  the  command  value  of 
400  g  is  shown  in  Fig.  6,  and  the  angle  of  attack  and  velocity  behavior  are 
shown  in  lig.  7.  The  angle-of-attack  growth  in  resonance  with  the  body-fixed 
trim  asymmetry  drives  the  vehicle  into  a  flat  spin  at  90-deg  angle  of  attack.  A 
significant  feature  of  the  control  system  response  is  the  relatively  low  fre¬ 
quency  of  the  c.g.  offset  or  moving-mass  oscillations  compared  with  the  char¬ 
acteristic  pitch  frequency  of  the  vehicle.  This  can  be  seen  from  Eq.  (7)  for 
the  linear  approximation  to  the  roll  angle  oscillations  relative  to  the  wind. 

For  small  oscillations  in  <j>,  Eq.  (7)  represents  a  harmonic  oscillator  with 
natural  frequency  ft  given  by 

ft  -  (nec)1/2  -  4(-£.)(i-)e]1/2  (32) 

St  X 

In  the  time  period  around  0.5  to  1  sec,  from  Figs.  5  and  7,  c  *  0.015 
in.,  0  »  12  deg  and  ui  «25  cps,  which  gives  for  ft  the  value  0.12  w  or  3  cps. 

This  agrees  well  with  the  oscillation  frequency  observed  in  Fig.  5.  The  c.g. 
offset  amplitude  required  for  control  is  approximately  ±  0.020  in. ,  from  Fig. 

5b.  This  requires  a  moving-mass  throw  weight  of  only  2.5  lb-ln.,  e.g.,  a  2.5- 
lb  mass  with  a  displacement  of  ±1  in.  or  a  1.25-lb  mass  with  a  displacement  of 
±2  in.,  for  the  122-lb  example  vehicle.  In  view  of  the  relatively  low  oscil¬ 
lation  frequency  required  for  this  mass,  the  energy  requirements  of  the  con¬ 
troller  are  minimal. 


Roll-rate  response,  c.g.  offset  movement,  and 
windward-meridian  angle 
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Fig.  7.  Recovery  angle  of  attack  and  velocity 
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V.  CONCLUSIONS 


The  concept  of  using  a  moving-mass  center~of-gravity  controller  to  limit 
roll  resonance  angle-of-attack  amplification  of  a  ballistic  reentry  vehicle 
has  been  demonstrated.  A  control  law  derived  from  a  first-order  linear 
approximation  to  the  nonlinear  rotational  equations  of  motion  near  resonance 
yields  a  stable  response  to  an  angle-of-attack  control  system.  The  first 
order  solution,  obtained  with  little  optimization,  is  an  indication  of  the 
stability  of  the  resonance  lockin  condition  driven  by  a  body-fixed  trim 
asymmetry  with  an  orthogonal  center-of-gravity  offset.  A  significant  feature 
of  the  resonant  motion  is  the  relatively  low  frequency  of  coupled  oscillations 
in  angle  of  attack,  roll  rate,  and  roll  angle,  about  the  steady-resonance 
condition.  This  low-frequency  motion  determines  the  response  requirements  of 
a  moving-mass  c.g.  controller,  and  the  energy  requirements  prove  to  be 
minimal.  The  resonance  control  system  has  been  demonstrated  with  an  applica¬ 
tion  to  ballistic  reentry  vehicle  recovery.  The  controller  can  limit  the 
angle  of  attack  and  normal  load  response  in  resonance,  while  providing  a  large 
angle-of -at tack-induced  drag  for  recovery.  A  digital-computer  simulation 
demonstrates  that  resonance  control  can  be  utilized  to  decelerate  a  ballistic 
reentry  vehicle  from  high  hypersonic  velocity  at  relatively  low  altitude  to 
subsonic  velocity  prior  to  impact.  The  vehicle  1b  ultimately  driven  into  a 
flat  spin  with  enormous  drag  deceleration,  while  limiting  the  normal  loads 
during  angle-of-attack  buildup  to  a  prescribed  value. 
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normal  acceleration 

normal  acceleration  command 

roots  of  Eq.  (15) 

radial  c.g.  offset 

axial  force  coefficient 

drag  coefficient 

lift  coefficient 

pitch  damping  derivative 
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normal  force  coefficient 

normal  force  derivative 

cN  qS/mu 
a 

aerodynamic  reference  (base)  diameter 
acceleration  due  to  gravity 
altitude 

altitude  at  recovery  initiation 
pitch  or  yaw  moment  of  inertia 
roll  moment  of  inertia 
feedback  gain,  (xaY.^/2 
feedback  gains 
vehicle  mass 
roll  rate 

roll  rate  at  recovery  initiation 
dynamic  pressure 
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earth  radius 

Laplace  transform  variable 
aerodynamic  reference  (base)  area 
time 

velocity 

velocity  at  recovery  initiation 
static  margin 
sin  ^  X 
path  angle 

path  angle  at  recovery  initiation 
roll  angle,  tt/2  -  4> 

CN  <S/I* 
a 

CN(0)qS/Ix 
angle  of  attack 

command  value  of  angle  of  attack 
pitch  rate  in  wind  coordinates 
tan  a /tan  0 

V1 

*•  * 

pitch  or  yaw  damping  parameter,  C  ■+* 

mq  a 

atmospheric  density 
cone  half  angle 

trim  angle  of  attack  asymmetry 
roll  angle  relative  to  wind 

meridian  angle  between  trim  asymmetry  and  c.g.  offset  (Fig.  1) 
precession  angle 
precession  rate 

undamped  natural  pitch  frequency 

1/2 

controller  frequency,  (n0c) 


30 


LABORATORY  OPERATIONS 


The  Laboratory  Operations  of  Tha  Aerospace  Corporation  la  conducting 
experimental  and  thaoratlcal  lnvaatlgatlona  nacaaaary  for  tha  avaluatlon  and 
application  of  aciantific  advancaa  to  new  military  concapta  and  systems.  Ver- 
aatlllty  and  flexibility  hava  baan  developed  to  a  high  dagraa  by  tha  laborato¬ 
ry  paraonnal  In  dealing  with  tha  nany  problaaa  encountered  in  tha  Nation' a 
rapidly  developing  apace  ayataaa.  Expertise  in  tha  la teat  aciantific  develop¬ 
ment!  la  vital  to  tha  accoapliahaent  of  taaka  related  to  theaa  problaaa.  The 
laboratorlaa  that  contribute  to  thle  raaaarch  arat 

Aarophyalca  Laboratory!  Aarodynaalca;  fluid  dynaalca;  plaaaadynaaica; 
chaalcal  klnatica;  engineering  aachanlca;  flight  d/naaica;  heat  tranafer; 
high-power  gaa  laaara,  contlnuoua  and  pulsed,  IR,  visible,  UVj  laser  phyalca; 
loser  raaonator  optica;  laser  affacta  and  countaraaaauraa. 

Chaalatry  and  Phyalca  Laboratory!  Ataoapharlc  reaction*  and  optical  hack- 
grounda;  radiative  tranafer  and  ataoapharlc  transmission;  thermal  and  atata- 
e pacific  reaction  rates  in  rocket  pluaas;  chaalcal  tharaodynaaica  and  propul- 
aion  chaalatry;  laaar  isotope  separation;  chaalatry  and  phyalca  of  pertidea; 
apaca  environmental  and  contamination  affacta  on  spacecraft  materials;  lubrlca 
tion;  surface  chaalatry  of  insulators  and  conductors;  cathode  materials;  sen¬ 
sor  materials  and  sensor  optica;  applied  laser  spectroscopy;  atomic  frequency 
standards;  pollution  and  toxic  materials  monitor log. 

Electronics  Raaaarch  Laboratory:  Electromagnetic  theory  and  propagation 
phaooeena;  microwave  and  semiconductor  device*  and  integrated  circuits;  quan¬ 
tum  (tlactronica,  laaara,  and  alactro-optlcs;  communication  adancas,  applied 
electronics ,  superconducting  and  electronic  device  phyalca i  mllllaatar-wave 
and  far~lnfrarad  technology. 

Materials  Sciences  Laboratory!  Development  of  new  material* j  coapoalta 
aatarlals;  graphite  and  caraalca;  polymeric  materials;  weapons  effects  and 
hardened  materials;  materials  for  electronic  devices;  dimens tonally  stable 
materials;  chemical  and  structural  analyses;  atraas  corrosion;  fatigue  of 
metals. 

Specs  8clancta  Laboratory!  Atmospheric  and  ionospheric  physics,  radia¬ 
tion  from  tha  atmosphere,  density  and  composition  of  tha  atmoaphara,  auroraa 
and  airgiow;.  magnatoepheric  physics,  cosmic  rsys,  generation  and  propagation 
of  plasma  waves  in  tha  nagnatoaphara;  solar  phyalca,  x-ray  astronomer;  tha  affacta 
of  nuclear  explosions,  magnetic  stoma,  and  tolar  activity  on  tha  earth's 
atmoaphara,  lonoaphera,  and  nagnetosphere;  the  affacta  of  optical,  electromag¬ 
netic,  and  particulate  radiations  in  apace  on  space  systems. 


